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We study thermalised configurations of SU(2) gauge fields by cooling. An analysis of the effect of cooling is 
presented and global and statistical information is extracted. 


1. Introduction 

It is highly desirable to have a qualitative un¬ 
derstanding of the mechanism underlying the oc¬ 
curence of non-perturbative phenomena like Con- 
Rnement, spontaneous chiral symetry breaking, 
the rj/ mass generation, etc., and a consequent 
estimation of the relevant physical quantities: 
string tension, topological susceptibility, < iptp >, 
etc. Many authors have attempted this program, 
giving a description of the physical Yang-Mills 
vacuum, and identifying what they believe are 
the relevant degrees of freedom to understand all 
these phenomena. The lattice formulation has 
been able of determining most of these physi¬ 
cal constants numerically out of the QCD la- 
grangian in fairly good agreement in most cases 
with the expectations. Hence, it seems that the 
lattice configurations which give rise to these re¬ 
sults contain the essential information to discrim¬ 
inate between the different mechanisms and mod¬ 
els. Some of these models are based in attribut¬ 
ing a fundamental role to a given class of classical 
configurations, and it should be possible a priori 
to identify them and check that they actually do 
as expected. 

In practice, the problem becomes very hard 
since lattice configurations are very noisy and it 
is hard to see any pattern in them. However, 
most authors agree that the noise is basically 
short ranged and does not affect long distance 
correlations, which lie in the heart of all non- 
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perturbative behaviour. The cooling method was 
introduced precisely as a method to get rid of the 
short range fluctuations while keeping the long 
range ones unaltered, and hence is the instrument 
used by most authors to attempt this program. 
Our work complies to this scheme. 

We present preliminary results based on the 
analysis of 56 configurations of SU(2) gauge fields. 
These configurations are generated by a heat bath 
method and Wilson’s action at /3 = 2.325 and a 
8^ X 64 lattice and with a m = (1, 1, 1) twist. 
Our configurations are separated by at least 7500 
Monte Carlo iterations and all our checks imply 
that they are fairly uncorrelated. The choice of 
parameters is based on our previous work [2] in 
which we studied SU(2) gauge fields for this twist 
value and with various /3 values and lattice sizes. 
There we saw that at the present values results 
have small finite size effects. 

2. The effect of cooling 

Since we are going to use cooling, our first con¬ 
cern is to know to what extent does cooling affect 
the results which are obtained. For that purpose 
we have carried on several checks and tests by 
applying the cooling procedure to known classi¬ 
cal configurations. We will actually be playing 
with various versions of cooling based of actions 
different from Wilson’s one. Our family of cool¬ 
ings [1] depends on one parameter e which is equal 
to 1 for Wilson’s action and to 0 for a classically 
improved action ( See Ref. 1 for details). Nega¬ 
tive values of e lead to opposite sign contributions 
to the action and are hence called over-improved. 
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This additional freedom is particularly important 
when the family of classical conhgurations con¬ 
tains zero modes, since then 0(a) corrections be¬ 
come leading. In particular, in Ref. 1 we saw 
that when considering sufhciently smooth lattice 
instantons, ordinary cooling causes the width to 
decrease, while for over-improved values the ten¬ 
dency is reversed. 

Our hrst test involved putting an isolated in- 
stanton on the lattice and studying the way in 
which it evolved with the different coolings. We 
have started with a lattice instanton extracted 
from our Monte Carlo data, and cooled it with 
€ = 1 until it became fairly small ( Thin Instan¬ 
ton). Then we reversed the sign of e and cooled 
for 1000 steps (Fat Instanton). Next, we cooled 
back with -fl until the instanton eventually dis¬ 
appeared. For each step we computed the quan¬ 
tity R = 67r^/£'™““^, where is the maxi¬ 

mum of the space integrated energy. For a con¬ 
tinuum instanton i? is equal to p. We also com¬ 
puted a series of improved versions of the topo¬ 
logical charge, which for smooth Reids reproduce 
the continuum formula to order The formula 
is based on computing the Chromo-Electric and 
Magnetic field strengths by a clover average of 
square plaquettes of sizes 1, 2, ...n. and then com¬ 
puting the integral of the scalar product in space 
and color of both fields: Q*-"). Fig. 1 displays the 
history profile, where the horizontal axis counts 
cooling steps and the kinks mark the points where 
e was reversed. Notice that for large enough in¬ 
stantons the width varies linearly. By varying the 
value of e we have found that the rate of growth 
or decrease depends linearly with e although with 
a different proportionality for positive and nega¬ 
tive e. 

For small instantons the behaviour is more com¬ 
plicated. We have analised the result of cooling 
with various e the Thin Instanton configuration. 
Actually all values larger that —0.3 lead to a de¬ 
struction of the configuration in a few tens of 
cooling steps. We conclude with the following 
panorama: Positive values of e lead to an even¬ 
tual destruction of all instantons in a calculable 
number of steps depending on the original size. 
For negative e values there is a minimal size be¬ 
low which they are unstable under cooling. Based 



Figure 1. R (thick line) and 4 ( top thin) 

and as a function of cooling step for a lattice 
instanton configuration. The value of e changed 
from 1 to -1 and back to 1. 


on these results we have chosen e = —0.3 as the 
value to be used in the remaining. 

It is also interesting to realize that due to the size 
variation, will also vary with cooling induc¬ 
ing a variation of the topological susceptibility. It 
is recommended to use which is very slightly 
p dependent and close to 1 even up to very small 
sizes. 


3. Global Quantities 

We have applied up to 50 cooling steps with e = 
—0.3 on the set of 56 configurations, and analysed 
the dependence of several observables with the 
number of cooling steps. We have measured the 
topological susceptibility associated to each of the 
clover averaged charges: 

Volume 

In Fig. 2 we show the results of as a func¬ 
tion of cooling steps, compared to the results ob¬ 
tained with Wilson’s action cooling. There is a 
sizable difference between both, which can be at¬ 
tributed to the destruction of instantons in the 
second case. Using instead of x*-^^ gives dif¬ 
ferences of the same magnitude. We suggest that 
the best choice would be to take e = —0.3 and cool 
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for a very large number of times. This should pro¬ 
duce a stable value of the susceptibility close to 
the one we found at 50 coolings. This value does 
not include the effect of instantons of sizes smaller 
than the Thin Instanton ( ~ 2 a), which have 
been eliminated by the cooling procedure. One 
should then estimate the contribution of these 
small instantons. In principle, asymptotic free¬ 
dom allows to determine the distribution of small 
enough instantons. 

Now we turn to other observables: The electric 
flux e sector energies, measured by looking at the 
correlation of spatial Polyakov loops of winding 
equal to e. We have measured the following quan¬ 
tities for all cooling steps up to 50 
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Figure 2. The topological susceptibility for 
e = —0.3 ( Riled squares) and e = 1 (filled cir¬ 
cles), and the quantities Ei(n() for rij = 3, 5 and 
7 (crosses, squares and tringles respectively) as a 
function of cooling steps. All quantities are given 
in appropiate fermi units. 


where e = |e| and Ce{nt) is the lattice correla¬ 
tion function at distance rij for Polyakov loops of 
flux e. The physical distance Ig = a{p)Ng is fixed 
by setting the infinite volume string tension to 
(T = bfm~^. In the infinite volume limit Yigirit) 
should tend to a for large rij independently of 
e. Our results are indeed quite compatible for 
all three fluxes, confirming again that finite size 
effects are small. Y^Prit) grows with rij for small 
values, but the data seems to flatten out for larger 
ones. Indeed, the result for rij = 5 is compatible 
with the larger distance ones, which have large 
errors. Taking the value of Se(5) as an estimate 
of the string tension, we could conclude that it 
decreases with cooling, being for 50 cooling steps 
at about 40 — 60% of its uncooled value. Since 
the configurations at fifty coolings are very closely 
given by a bunch of self-dual (instantonic) lumps, 
it seems that these objects account for at least 
50% of the value of the string tension. Of course 
these instantonic structures must be highly cor¬ 
related in order to give rise to a finite string ten¬ 
sion. In Ref. 2 we gave an scenario in which this 
can occur. Some work is needed to clarify the 
reason for the slow decrease that is still observed 
as a function of the number of coolings. 


4. Statistical Distributions 

We have performed a number of statistical 
analysis on our data sample, which cannot be 
presented in such short space. We found that 
the electric and magnetic field densities are highly 
correlated, but there is a slight excess of magnetic 
energy. Identifying energy peaks with instantons, 
we obtained that the distribution of widths is 
fairly fiat in contrast with the results of Ref. 3. 
The distribution in space is also fiat beyond 8 a, 
showing both attractive and repulsive region be¬ 
low. 
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